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We introduce a grating assisted tunneling scheme for tunable synthetic magnetic fields in pho¬ 
tonic lattices, which can be implemented at optical frequencies in optically induced one- and two- 
dimensional dielectric photonic lattices. We demonstrate a conical diffraction pattern in particular 
realization of these lattices which possess Dirac points in fc-space, as a signature of the synthetic 
magnetic fields. The two-dimensional photonic lattice with grating assisted tunneling constitutes 
the realization of the Harper-Hofstadter Hamiltonian. 

PACS numbers: 42.50.Xa, 42.82.Et, 03.65.Pm 



Synthetic magnetism for photons is a unique tool for 
the manipulation and control of light, and for the de- 
sign of novel topological phases and states in photon- 
. Topological photonics is a rapidly growing 
advancing in parallel to analogous efforts in ul¬ 
tracold atomic gases [lallll: inspired by the development 
of topological insulators in condensed matter physics [l^ . 
One motivating aspect of topological photonic systems 
is the existence of unidirectional backscattering immune 
states n , which are robust to imperfections, and thus 
may serve as novel waveguides and for building integrated 
photonic devices. The first experimental observations 
of such edge states were in the microwave domain, in 
magneto-opticalphotonic crystals 0, theoretically pro¬ 
posed in Refs. [1H3- In the optical domain, imaging of 
topological edge states was reported in Floquet topolog¬ 
ical insulators, inmlemented in modulated honeycomb 
photonic lattices (^, and in the two-dimensional array 
of coupled optical-ring resonators The strategies 
for obtaining synthetic magnetic/gauge fields and topo¬ 
logical phases for optical photons are closely related to 
the system at hand. In systems of coupled optical res¬ 
onators, the strategy is to tune the phase of the tun¬ 
neling between coupled cavities for example, by 

using link resonators of different length in or time- 
modulation of the coupling [^. Photonic topological 
insulators were also proposed in superlattices of meta¬ 
materials with strong magneto-electric coupling [l^ . In 
2D photonic lattices (waveguide arrays), pseudomagnetic 
fields have been demonstrated by inducing ’strain’ in op¬ 
tical graphene 0- By modulating ID photonic lattices 
along the propagation axis, one can choose the sign of 
the hopping parameter between neighboring sites m, 
whereas topological states were achieved using ID pho¬ 
tonic quasicrystals 12|. However, even though photonic 
lattices possess a great potential for exploring synthetic 


magnetism and topological effects, a viable scheme for 
arbitrary design of the phases of the complex tunnel¬ 
ing matrix elements still needs to be developed. Here 
we introduce one such scheme termed grating assisted 
tunneling, and propose its irnplementation in optically 
induced photonic lattices 119142311 We demonstrate the 
conical diffraction pattern |24l . l25l| in particular realiza¬ 
tions of one- (ID) and two-dimensional (2D) square pho¬ 
tonic lattices with grating assisted tunneling. The latter 
constitutes demonstration of the the (’single particle’) 
Harper-Hofstadter Hamiltonian (HHH) [2^, in these 
systems. 

The proposed method is inspired by the so called 
laser assisted tunneling scheme, which was implemented 


in optical lattices with ultracold atoms [28 



scheme is based on theoretical proposals in Refs, 
subse quen tly developed to experimentally realize the 
HHH [29|, Isdl . a nd staggered magnetic fields in optical 
superlattices [28|. Many of the results obtained with ul¬ 
tracold atoms are viable in photonic lattices, and per¬ 
haps even more of them are feasible, because heating by 
spontaneous emission that is present in ultracold atomic 
systems, is absent in photonic systems. 

We consider the paraxial propagation of light in a 
photonic lattice defined by the index of refraction n = 
no + 6n{x, y, z) {Sn <C no), where no is the constant back¬ 
ground index of refraction, and Sn(x, y, z) describes small 
spatial variations, which are slow along the propagation 
z-axis. The slowly varying amplitude of the electric field 
ipix, y, z) follows the (continuous) Schrddinger equation: 
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here, = d'^ jdx^ jdy^, and k = 27rno/A, where A is 
the wavelength in vacuum. From now on we will refer to 
Sn(x, y, z) as the potential or index of refraction. In the 
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simulations we use ng = 2.3, corresponding to the sys¬ 
tems that were used to implement the optical induction 
technique [2^[2l|, and A = 500 nm. 

For clarity, let us first present the grating assisted tun¬ 
neling method in a ID photonic lattice. If the potential 
is a periodic lattice, 6nL{x) = SnL{x + a), where a is the 
lattice constant, and if the lattice is sufficiently ’deep’, 
the propagation of light can be approximated by using a 


discrete Schrodinger equation (e.g., see 
atoms): 


for ultracold 


discrete Schrodinger equation [23 
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where Jx quantifies the tunneling between adjacent 
waveguides and is the amplitude at the mth lat¬ 

tice site, i.e., waveguide. To create a synthetic magnetic 
field, we adopt the strategy to tune the phase of the 
tunneling between lattice sites. For this ID lattice, we 
seek a scheme which effectively renormalizes Jx to get 
Kxeyzp(i4>m), where (j)m denotes the phase for tunneling 
from site m to site m-|-l. The scheme is illustrated in Fig. 
[TJ In Fig. [ija) we show the potential 5n{x, z), which can 
be modeled as a discrete lattice with complex tunneling 
parameters Kx shown in Fig. [IJb). 

Let us gradually explain the idea behind the variation 
of the index of refraction as in Fig. (TJa). In Fig. [IJe), 
we show the propagation of intensity in a continuous ID 
model in the lattice potential Sn^^x) = Sn^Q cos^{'kx/ a), 
with Stilo = 4 x 10“^, a = 

10 /rm. We see the usual diffraction pattern for a spatially 
broad excitation covering several lattice sites 23|. Next, 
suppose that we introduce a linear gradient of index of re¬ 
fraction along the x direction Sririx) = —rjx in addition 
to the lattice potential, such that Sn = SnL(x) + 5nT(x). 
For a sufficiently large tilt, the tunneling is suppressed. 
This can be seen from Fig. [Ijf) which shows the propaga¬ 
tion of intensity in a tilted potential with rj = O.lSriLQ/a; 
the tilt should be smaller than the gap between the first 
two bands. Finally, let us introduce an additional small 
grating potential at a small angle 9 with respect to the 
z axis, Sncixjz) = Sncocos^{{qxx — Kz)j2)^ such that 
5n(x,z) = 5nLix)+SnT{x)+6noix, z). This total poten¬ 
tial Sn{x, z) is illustrated in Fig. [T])a). The ’frequency’ k 
is determined by the angle 9 of the grating with respect to 
the z-axis, which is chosen such that k = pak/uQ and the 
grating forms a z-dependent perturbation resonant with 
the index offset between neighboring lattice sites pa (Fig. 
[ija)). The grating restores the tunneling along the x- 
axis, hence the term grating assisted tunneling. Restored 
tunneling is seen in Fig. [ijg) which shows diffraction for 
identical initial conditions as in Figs. [Ije) and (f); the 
grating parameters are Qx = tt/o and Snco = Q.ISulq. 

However, the diffraction pattern is drastically changed. 
In order to interpret it, we point out that, for resonant 
tunneling where k = pak/no and a sufficiently large tilt 
(J p), z-averaging over the rapidly oscillating terms 
shows that the system can be modeled by an effective 


. diprr. 

dz 


= -{Kxe^<^^-^^m-i + ( 3 ) 


where 4>m = q ■ Rm = qxxna. In Fig. [Ijg) we used 
qx = 7r/a, i.e., (j)m = rmr. Such a discrete lattice is il¬ 
lustrated in Fig. [ijc); its dispersion having a ID Dirac 
cone a,t kx = 0. For a wavepacket that initially excites 
modes close to kx = 0, the diffraction in the discrete 
model ©yields the so-called (ID) conical diffraction 
pattern [llil24|, as illustrated in Fig. [TJh). The initial 
conditions for propagation in the discrete model corre¬ 
sponds to the initial conditions in the continuous system, 
V’m(O) = \fTe~^™‘l'^'^ , and Kx = 0.053 mm“^. Thus, we 
interpret the diffraction pattern in Fig. [TJg) as ID coni¬ 
cal diffraction, a signature of the discrete model depicted 
in Fig. (He). A comparison of the discrete [Fig. [TJh)] 
and the realistic continuous model [Fig. HKg)], clearly 
shows that we can use grating assisted tunneling to tune 
the phases of the tunneling parameters in the discrete 
Schrodinger equation, thereby realizing synthetic mag¬ 
netic fields. 

Before proceeding to 2D systems, we discuss the ampli¬ 
tude of the tunneling matrix elements Kx as a function of 
the strength of the grating dneo- Figure [T][d) shows Kx 
(green squares) and Jx (blue circles) versus Sneo, where 
Jx corresponds to the potential which includes the lat¬ 
tice and the grating, but no tilt. The amplitudes Jx and 
Kx are obtained by comparing the diffraction pattern of 
the discrete with the continuous model, and adjusting Jx 
and Kx until the two patterns coincide, as in Figs. [TJg) 
and (h). Our results in Fig. [TJd) are in agreement with 
those in ultracold atoms [e.g., see Fig. 3(a) in Ref. [^]. 

The extension of the scheme to 2D lattices is 
straightforward. We consider a square lattice, Sul = 
6nLo{cos^{'!Tx/a) + cos"^{ny/a)) , the tilt in the x di¬ 
rection, Snx^x) = —px, and the grating which has 
the form 5nc{x^y, z) = Sneo cos^{{qxX + qyp — kz)I2). 
Propagation of light in the total potential 5n{x, y, z) = 
Sni^{x, y) + SnT{x) + 5nc{x, y, z) can be modeled by the 
discrete Schrodinger equation (the derivation is equiva¬ 
lent to that in Ref. [1^ for ultracold atoms): 


. dlprn.r. 


l—r— = -{Kxe^’>’^-^-^lljm-l,n + Kxe 
dz 


+ Jylj} 

m,n— 1 T 'Ty'(/^m,n-t-l); 


( 4 ) 


where 4>m,n = q • R-m.n = qxXna + qyua. Note that the 
tunneling along y does not yield a phase because there is 
no tilt in the y direction; the tunneling amplitude along 
y depends on the depth of the grating, as illustrated in 
Fig. [TJd) with blue circles. 

In order to demonstrate that the propagation of light 
in the continuous 2D potential Sn{x, y, z) is indeed equiv¬ 
alent to the dynamics of Eq. a, we compare prop¬ 
agation in the discrete model a, with that of the 
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FIG. 1. (color online) Illustration of the grating assisted tunneling scheme, (a) Sketch of the spatial dependence of the index of 
refraction Srihix, z) = Sul + Sut + Sna, which creates a synthetic magnetic field in a photonic lattice. A ID photonic lattice 
[Sul = SriLo cos^{nx /a)] is superimposed with a linear gradient index in the x direction (Sut = —rix, orange dotted line), and 
an additional small grating potential [Sna = Snao cos^{{qxX — kz)/2)], at a small angle {6 is on the order of 1°) with respect 
to the a axis. In the xz plane we plot the projection of Stil (blue bold lines), and the grating (green dashed lines), (b) A 
discrete lattice with complex tunneling matrix elements between sites, , and spatially dependent phases (j>m, can model 

the system in (a), (c) A particular choice of 6 (see text), yields a lattice with (f>m = nm, and a dispersion with a ID Dirac 
cone at kx = 0. (d) The amplitude of the tunneling Kx as a function of the strength of the grating Snao (green squares). The 
tunneling amplitude Jx versus Snao for a system without the tilt. (e,f,g) Numerical simulation of the evolution of a wavepacket 
that initially excites modes close to = 0 in the continuous ID photonic lattice, (e) Diffraction in a periodic photonic lattice 
SriL- (f) Propagation in the tilted system, Sul + Sut, shows that the tunneling (diffraction) is suppressed, (g) The tunneling 
is restored by an additional grating potential with qx = ir ja. (h) Propagation of the wavepacket in the corresponding discrete 
model illustrated in (c). 


continuous equation O- The lattice parameters are tion for 2 ; = 162 mm. The two concentric rings are a clear 


SriLQ = 4 X 10 
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= 13 /rm; the tilt is given by 


r] = O.lSriLo/a; the grating is defined by Snco = O-lTdriLo 
and Qx = —Qy = 7r/a, which yields (prn,n = (w ~ 
and K is chosen to yield resonant tunneling. The dis¬ 
crete lattice which corresponds to this choice of phases 
is illustrated in Figure [UJa). It has two bands, /3 = 

Kx sin^(fca;a) -h Jy cos^(fcya) (/3 is the propagation 
constant), touching at two 2D Dirac points at (fc^,, ky) = 
(0, ±7r/2a) in the Brillouin zone 3^, as depicted in 
Fig. [2Db). Suppose that the incoming beam at 2 : = 0, 
excites the modes which are in the vicinity of these two 
Dirac points. Around these points, for a given k = k/fc, 
the group velocity, Vk/3(fca;,/cj,), is constant. Thus, the 
beam will undergo conical diffraction, which has been 
thoroughly addressed with Dirac points in honeycomb 
optical lattices [ 2 ^. To demonstrate this effect, we con¬ 
sider the propagation of a beam with the initial pro¬ 
file given by ip{x,y,0) = •\/Tcos(2/7r/2a) exp(—(a;/3a)^ — 
{y/Sa)"^) in the 2D potential 6n{x,y,z) (the cosine term 
ensures that we excite modes close to the two Dirac 
points). In Fig. [2Kc) we show this beam after propaga- 


evidence of conical diffraction [2J, [2^ . We also compare 
this with the propagation in the discrete model (jT]), with 
tunneling phases plotted in Fig. [2Ia); = 0.11 mm“^, 

Jy = 0.14 mm“^. The results of the propagation in the 
discrete model are shown in Fig. [5Dd). It is evident that 
for our choice of the tilt and the grating, we have effec¬ 
tively realized the lattice plotted in[2Ka). This is in fact 
a realization of the HHH for a = 1/2, where a is the flux 
per plaquette in units of the flux quantum 2^, ^ . 


For the experimental implementation of the scheme, 
we propose the so-called optical induction technique in 
photosensitive materials, which can be implemented in 
photorefractives 1^23|. In these systems, both the lat¬ 
tice 5nL{x,y) and the grating potential Sncix^y, z) can 
be obtained in a straightforward fashion by using inter¬ 
ference of plane waves in the medium (2^. The lattice 
constant a, the grating parameter q, and hence the hop¬ 
ping phase (prn,n = q ■ are tunable by changing 

the angle between the interfering beams. This could en¬ 
able manipulation of the phases of the tunneling matrix 
elements in real time I20II. 
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FIG. 2. (color online) Grating assisted tunneling and coni¬ 
cal diffraction in a square photonic lattice, (a) Sketch of the 
2D lattice with grating assisted tunneling along the x direc¬ 
tion. The resulting nontrivial hopping phases tt and 0 are de¬ 
noted with dashed and solid lines, respectively. A wavepacket 
that makes one loop around the plaquette accumulates the 
phase TT. (b) The lattice possesses two 2D Dirac cones at 
{kx,ky) = (0, ±7r/2a) in the dispersion l3{kx,ky), where f3 is 
the propagation constant, (c) Intensity of a beam, which ini¬ 
tially excites modes in the vicinity of the Dirac points, after 
propagation for z = 162 mm. The intensity has two con¬ 
centric rings corresponding to the conical diffraction pattern, 
(d) Simulation of the diffraction pattern in the discrete model 
corresponding to the lattice in (a), which also exhibits conical 
diffraction (see text for details). 


The most challenging part of the implementation ap¬ 
pears to be creation of the linear tilt potential. To ob¬ 
serve the conical diffraction effects, one needs a linear 
gradient over the ^ 20 lattice sites, which implies a to¬ 
tal index tilt An = 20T]a ~ 20 x O.lJriLo ~ 8 x 10“'^. 
In principle, the linear tilt could be achieved by using a 
spatial light modulator. Another possibility which could 
create a linear tilt is to use crystals with linear depen¬ 
dence of the index of refraction on temperature, and a 
temperature gradient across the crystal. 

However, it should be emphasized that if one uses a su¬ 
perlattice for 6nL, rather than the linear tilt potential, as 
in Ref. [2^ for ultracold atomic gases, the grating assisted 
tunneling scheme proposed here would straightforwardly 
yield staggered synthetic magnetic fields for photons. The 
achievement of such staggered fields is straightforward 
with optically induced lattices proposed here. 

Before closing, we emphasize that the spatial phase be¬ 
tween the periodic lattice potential Sul and the grating 
dnc is a relevant parameter. This point has recently been 
examined in detail in the context of ultracold atoms ISSlI . 


To translate it to our system, suppose that in our ID sys¬ 
tem, we shift the grating along the a:-direction, such that 
Snc = dnco cos^{{q^x — nz + ^)/‘2). This simply shifts 
the phases in the discrete lattice [shown in Fig.[T][c)] such 
that (j)m is replaced by + which moves the ID Dirac 
point in fc^-space by Akx = ^/a. We have numerically 
verified that this indeed happens. 

In conclusion, we have introduced a grating assisted 
tunneling scheme for tunable synthetic magnetic fields 
in photonic lattices, and proposed its implementation 
at optical frequencies in optically induced one- and two- 
dimensional dielectric photonic lattices. As a signature 
of the synthetic magnetic fields, we have demonstrated 
the conical diffraction pattern in particular realization of 
these lattices [shown in Fig.[T][c) and Fig.|2l[a)]. The two- 
dimensional photonic lattice with grating assisted tunnel¬ 
ing constitutes the realization of the HHH. The scheme is 
well suited for realization of staggered synthetic magnetic 
fields in photorefracives. We envision that this proposal 
will open the way to other studies of light propagation 
in photonic lattices with complex tunneling matrix ele¬ 
ments, including nonlinear propagation, solitons and in¬ 
stabilities in synthetic magnetic fields, and creation of 
synthetic dimensions in photonic lattices by employ¬ 
ing synthetic fields. 
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Note added. A few days before this paper was sub¬ 
mitted to arXiv, a paper by a different group [s^ ap¬ 
peared on the arXiv. It is also inspired by laser as¬ 
sisted tunneling in cold gases, and experimentally demon¬ 
strated modulation-assisted tunneling in laser-fabricated 
photonic Wannier-Stark ladders. 
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